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ABSTRACT 

We have investigated how envelope poUution by icy planetesimals affects the critical 
core mass for gas giant formation and the gas accretion time-scales. In the core- 
accretion model, runaway gas accretion is triggered after a core reaches a critical core 
mass. All the previous studies on the core-accretion model assumed that the enve- 
lope has the solar composition uniformly. In fact, the envelope is likely polluted by 
evaporated materials of icy planetesimals because icy planetesimals going through 
the envelope experience mass loss via strong ablation and most of their masses are 
deposited in the deep envelope. In this paper, we have demonstrated that envelope 
pollution in general lowers the critical core masses and hastens gas accretion on to 
the protoplanet because of the increase in the molecular weight and reduction of adia- 
batic temperature gradient. Widely- and highly-polluted envelopes allow smaller cores 
to form massive envelopes before disc dissipation. Our results suggest that envelope 
pollution in the course of planetary accretion has the potential to trigger gas giant 
formation with small cores. We propose that it is necessary to take into account enve- 
lope pollution by icy planetesimals when we discuss gas giant formation based on the 
core-accretion model. 
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1 INTRODUCTION 

More than 540 exoplanets have been discovered so fai0. Ded- 
icated planet surveys have revealed the diversity of giant 
planets outside the Solar System. Many Sun-like stars are 
known to harbor close-in giant planets called hot-Jupiters, 
in contrast to our Sun. Recently, direct imaging has shed 
light on the existence of distant extrasolar giant planets 
such as HR8799b, c. d, e (|Marois et al.|[20oi.l2010l'). Fomal- 
haut b (|Kalas et al.ll200^ ). and Beta Pic b ( Lagrange et al.l 
I2OO9I ). A large number of exoplanets enable us to consider 
planet formation around various types of star through a 
statistical approach, which is often called population syn- 
thesis. In the population synthesis, planetary systems are 
built by putting together pieces of formation processes (e.g., 
solid and gas accretion, planetary migration and so on) in 
a Monte-Carlo way. Nowadays several groups h ave worked 
on the population synthesis independently (e.g. Ida fc Lin 
200i; lAlibert. Mordasini. fc Benjboill: [Kennedy fc Kenvon 



indicates that our basic picture of planetary formation is 
supported by observed exoplanets. On the other hand, it 
is also true that different population syntheses based on 
different assumptions as to planetary accretion and disc 
properties reproduce the same observed distribution. This 
highlights the need to understand underlying formation pro- 
cesses of giant planets more exactly. 

The core-accretion model is one of the promising 
models for giant plan et for mation (e.g. Mizund Il980l : 



iBodenheimer fc Pollack Il986l : IPoUack et al.l 



1996l h In this 



2008; M iguel. Guilera. fc Bruninill201ll '). Those population- 
synthesis studies succeed in reproducing the planetary mass- 
period distribution of detected extrasolar giant planets. This 
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model, once a core reaches a critical core mass through 
planetesimal accretion, runaway gas accretion is triggered 
and the envelope mass increases rapidly. The typical value 
of the critical core mass is thought of as being > lOMm 
JPoUack et al.l Il996l : iFortier. Benvenuto. fc Bruninil I200I 
l2009t ) to form a massive env elope within the disc l i fetime 
suggested by o b serva tions (|Haisch. Lada. fc Ladal 1200 ll : 
iHernandez et all |2007| ) . This is, however, incompatible 
with the inferred mass of Jupiter's core of < IQM^ 
(Isaumon fc Guillot|[2003 ). The critical core mass depends 
on plan etesimal accre ti on rate and opacities in the en- 
velope ( Mizund Il98d : llkoma. Nakazawa. fc Emoril l200d : 
iRafikovl 2006h . Fast accretion of planetesimals leads to 
form gas giants with large cores. Although slow rates of 
planetesimal accretion lead to smaller critical cores of 
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< IOM0, gas accretion on to small cores tends to be 
long. Reduction of grain opacities in the envelope can not 
only lower critical core masses but also hasten gas accre- 
tion on to the protoplanet jlkoma. Nakazawa. fc Emoril 
2OO0I: Hubickvi. Bodenheimer. fc Lissauer] 20051 : 



Papaloizou fc NelsonI l2005l: i Hori fc Ikomal l2010h . For 
example, Hori fc Ikomal ( 201C ) showed that core masses 



needed to form massive envelopes before disc dissipation are 
as small as IMq in the extreme case of grain-free envelopes. 
In fact, recent calculations on the dynamical behaviour of 
dust grains in the accreting envelope demonstrate that grain 
opacities in the envelope can be reduced to be on average 
1% o r less of those in the protoplanetary disc (p^odolakl 



1/0 or less ot tnose m tne protopianetary gisc [l-roaoiaj' 
l2003l :[ Movshovitz fc Podolakll200i : ]Movshovitz et al.ll2010^ 



The reason is that small grains initially suspended in the 
outer envelope quickly grow large in size and then settle 
down into the deep envelope where temperature is high 
enough so that grains evaporate. Their results support that 
reduction of grain opacities can be a feasible idea that 
relaxes the problem of the formation of gas giants with 
small cores. 

Another factor that affects the critical core mass is 
considered in this paper. All the previous studies assumed 
that the composition of gas is solar throughout the en- 
velope. In fact, the envelope is likely to be polluted by 
evaporated materials of icy planetesimals. Once a proto- 
planet grows up to be Mars-size, it becomes difficult for 
the protoplanet to capture planetesimals because the pro- 
toplanet enhances their random velocities through gravi- 
tational scatterings by itself. On the other hand, disrup- 
tive collisions between planetesimals with high relative ve- 
locities take place frequently and yield a large number 
of 10 to 100 m-sized fra gments unless the planetesimals 
are larger than 100 k m (|lnaba. Wetherill. fc Ikomal l2003l : 
iKobavashi et all bOld V The protoplanet can capture such 
small-sized fragments because of the enhanced collisional 
cross-section b etween the envelope and those fragments 
due t o gas drag (|lnaba fc Ikomall2003l : lBenvenuto fc Bruninil 
I2OO8I ). The fragments going through the envelope experience 
both melting and evaporation and lose most o f their masses 
in th e deep envelope before reaching the c ore (| Pollack et al.l 
ll986l : |Podolak. PoUack. fc Revnoldslll988l ). As a result, en- 
velope pollution by evaporated materials of icy planetesi- 
mals occurs inevitably in the course of planetary accretion. 
In addition, evaporation (or erosion) of the core itsel f is als o 
expected to occur, as pointed out bv lLissauer et all (|l995h . 
This fact motivates us to consider giant planet formation 
taking into account the effects of the heavy-element enrich- 
ment in the envelope. 

Approximate solutions for the critical core mass analyt- 
ically derived by previous studies tell that envelope pollution 
has significant impacts on the critical core mass: 

(a) Wholly-radiative envelopes (|Stevensonlll982l ): 

Merit oc fi- ^kHiJ (In R^ut)'^, (1) 

(b) Wholly convective envelopes (jWuchter ill 19931 ): 



Merit oc ^ 2 r2 



i \/ri - 4/3 



Pou 



(2) 



(Pi -1)2 

where Merit is the critical core mass, the molecular weight 
K the Rosseland mean opacity. Me the planetesimal accre- 



tion rate, and Fi the adiabatic exponent. The subscript 
"out" means the value at the outer edge of the envelope. 
The envelope pollution causes increase in ^ and reduction of 
Vad (or reduction of Fi). The reduction of Vad comes from 
dissociations of several kinds of molecules. The two formu- 
lae indicate that both effects reduce the critical core mass. 
This is because those effects lead to reduction of the local 
pressure gradient that supports the gravity. On the other 
hand, the envelope pollution also causes rise in gas opacities 
as a negative effect. The rise in gas opacities, mainly II2O 
that absorbs infrared radiation, increases the critical core 
mass by increasing the temperature gradient that results in 
an increase in the pressure gradient. 

In this study, we show that envelope pollution by icy 
planetesimals lowers critical core masses in most cases and 
hastens gas accretion on to the protoplanet despite the rise 
in gas opacities. We present numerical procedures in Sec- 
tion 2, including our modelling for envelope pollution by icy 
planetesimals. Results of critical core masses and gas ac- 
cretion time-scales are shown in Section 3 and Section 4, 
respectively. Discussions and conclusion are in the last two 
sections. 



2 NUMERICAL PROCEDURE 

2.1 Structure and Evolution of the Envelope 

We suppose a spherically-symmetric protoplanet that con- 
sists of a gaseous envelope and a rigid core with a constant 
density of 3.2 g cm~^. The envelope is assumed to be in hy- 
drostatic equilibrium. Its structure is simulated with a set 
of the following basic equations (e.g. lKippenhahn fc Weigerj 
that includes the equation of hydrostatic equilibrium, 

dMr 47rr4 ' ^ ' 

the equation of mass conservation, 
dr _ 1 

the equation of heat transfer, 

dT T dP „ 



(4) 



dMr P dMr ' 

and the equation of energy conservation. 



dLr 



T 



dS 



(5) 



(6) 



dMr dt 
where Mr is the mass interior to radius r, P the pressure, 
T the temperature, p the density, S the specific entropy, 
Lr the energy fiux crossing at a sphere of radius r, V = 
dlogT/dlog P, G the gravitational constant, eaec the specific 
energy generation by planetesimal accretion, and t the time. 

The dominant mechanism of heat transfer is chosen by 
the Schwartzschild criterion for convective instability: 



dluT 
dlnP 



Vad < Vr 



(7) 



where Vad is the adiabatic temperature gradient and Vrad 
the radiative temperature gradient given by 

3 uPLr 



Vrad ~ 



l&iracG MrT^ ' 



(8) 
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where c is the velocity of light, a the radiation density con- 
stant, and K the Rosseland mean opacity. Opacities both of 
dust grains and gas and the equation of state for gaseous 
components that we use are mentioned in Section [2.21 

The boundary conditions axe given at the core surface 
and the outer edge of the envelope as follows: 



3M„ 



47rpcorc 



1/3 



and 1/ = I/c 



at Mr = Mc, 



and 



T — Tdisc and p = pdisc at Mr 



(9) 



(10) 



where Mcorc and pcoro are the core mass and density, re- 
spectively, I/corc the luminosity at the core surface, Mp the 
protoplanet's total mass, and pdisc and Tdisc are the density 
and temperature of the disc gas, respectively. The total mass 
of the planet, Mp, is the mass of gas contained inside the 
outer edge of the envelope, which is defined by the smaller 
of the Bondi radius (-Rb) and the Hill radius (i?Hiii): 



R GMp 

tlB ~ 2 — ' ^Hill ~ dp 



Afp 



1 1/3 



3(A/p 



(11) 



where Cs is the sound velocity, ftp the semimajor axis of 
the protoplanet, and Af, the mass of the central star. In 
this paper, we use ap = 5.2AU, M, — IMq, Tdisc = 
150K , and Pdisc ~ 5.0 x 10~^^ g cm~"^, unless otherwise 
noted. 

We assume that the energy release by planetesimals oc- 
curs in a narrow region on top of the core (i.e., face = 0). 
This assumption seems to be incompatible with the situa- 
tion that planetesimals evaporate and deposit their kinetic 
energy before reaching the core. However, such simplifica- 
tion affects the results little, because the energy deposition 
of planetesimals occurs mostly in the deep, convective enve- 
lope, the structure of which is insensitive to where energy 
deposition occurs inside itself. On the other hand, we neglect 
the effect of sinking of dissolved materials of icy planetesi- 
mals in this study. Gravitational potential energy released 
by sinking of ablated materials makes a contribution to slow 
down gas accretion on to a proto-gas giant as pointed out by 
iPoUack et all (| 19961 '). Whether the energy release by sinking 
of ablated materials occurs effectively in a wide region where 
planetesimals slow down via ablation is important for forma- 
tion time-scales of gas giants. It will be considered in future 
calculations based on trajectories of planetesimals through 
a protoplanet's envelope. 

2.2 Envelope Pollution by Icy Planetesimals 

In this study, we consider that the envelope consists of 
two components, materials of disc gas (component A) and 
those of icy planetesimals (component B). The chemical 
compositions of component A and B are assumed to be 
solar and comet-Halley-like, respectively. The mass frac- 
tions of hydrogen, helium, and the others in the solar abun- 
dances are Xa = 0.711, Ya = 0.274, an d Za = 0.015, 
respectively (|Lodders. Flame, fc Gailll2009l ). while those in 
comet Halley are Xb = 0.06, Yb ^ 0.00, and Zb = 0.94 
(|Mumma. Weissman. fc Stern|[l993 '). 

We consider that the envelope has two-layer structure 
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Figure 1. Schematic picture of the envelope polluted by icy plan- 
etesimals. The envelope has two-layer structure: the non-polluted 
(Za) and the polluted (Zh) layers, where Za and Z]^ are the mass 
fractions of heavy elements for the disc gas and the mixture of 
the disc gas and icy planetesimals, respectively. The boundary be- 
tween the upper and the lower layers is defined by the homopause 
temperature (Tj,). 



(see Figure [T]). The upper layer consists only of component 
A because disc gas flows into the outer part of the enve- 
lope. The lower layer consists of component B in addition 
to A. The reason is as follows. While going through the en- 
velope, planetesimals experience mass loss via strong abla- 
tion. Most of their m asses arc deposited in the deep envelope 
(|Follack et al.lll986l : [Fodolak, Follack, fc Rcvnolds 198^ ) . In 
the lower envelope, convection occurs because of molecu- 
lar dissociation and high gas opacity due to the bound-free 
absorption of hydrogen. The evaporated materials are thus 
stirred effectively by convection, so that the chemical com- 
position becomes uniform in the lower layer of the envelope. 
Such a boundary between convective and radiative regions 
(i.e., "tropopause") corresponds typically to a temperature 
of 2000K. 

However, the boundary is likely to be above the 
tropopause. Even in radiative regions, eddy diffusion car- 
ries the evaporated materials of icy planetesimals upward 
like in the homosphere extending above the tropopa use on 
the present Earth (e.g. IChamberlain fc Huntenlll987l ). Also, 
pollution occurs even in the upper envelope because the tem- 
perature there is enough to evaporate ice. While only a small 
fraction o f the mass of planetesimals evaporate s in the upper 
envelope (|Fodolak. Follack. fc Revnoldsll 19881 ) , the envelope 
is itself tenuous in such a region, so that the upper envelope 
can be also enriched in heavy elements. Thus, we regard 
the boundary between the two layers as a free parameter 
in this study. To do so, we introduce "homopause temper- 
ature", Th; the values of Th that we use in this study are 
300, 500, 1000, and 2000K. 

Thermodynamic quantities of the envelope gas are cal- 
culated under the assumption of chemical equilibrium. Al- 
though chemical equilibrium may not be achieved in the 
region of low temperature because chemical reactions pro- 
ceed slowly, this problem is beyond the scope of this pa- 
per. The following 13 constituents are considered, namely, 
H, He, C, O, Ha, O2, CO, CO2, H2O, CH4, H+, O", and 
e~ . Thermodynamic quantities are calculated from the ther- 
modynamic potentials. The relevant physical quantities are 
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given in NIST-JANAF the rmochemical tables (|Chaselll99i : 
lOtt fc Boerio-Goatesir2000l ). 

We determine the mass fractions of carbon (Z^) and 
oxygen (Z^) in the lower envelope in such a manner that 
the sum of the two fractions is equal to Zh and the ratio is 
conserved: 



Zh + Zh , 



7° 



e)Z^ + eZ^, 



(12) 



e)z2 



(13) 



where e is the mixing ratio of component B. The mass frac- 
tions of hydrogen and helium in the mixture are described 
by e as well: 

Yh = {l^e)YA + eYB. 

The ratio of to is 0.004 : 0.011, while that of Z§ to 
Zg is 0.24 : 0.69. In this study, we deal with Zh, which is 
equivalent to e, as a free parameter for simplicity because 
the actual value of Zh depends on various processes such as 
ablation efficiency of icy planetesimals and their entry ve- 
locities and should be determined in a complicated manner. 

The gas opacity is derived from opacity tables with 
different values of X and Z given by J. Ferguson 
l|Alexander fc FergusonI 1 19941 ). Their calculations include 
other heavy elements besides H, He, C, and O. All the 
elements except H and He are assumed to be in the solar 
abundances. This is reasonable because the abundances of 
heavy elements in com et Halley are known to be similar t o 
the solar abundance s (iMumma. Weissman. fc Sternlll993l ). 
The opacity tables of [Alexander fc FergusoiT i 19941 ) are pro- 



vided as a function of logQ = p/Tg , where p is the density 
in g cm~^ and Te the temperature in million K. If log Q 
exceeds the available range of the opacity tables, namely, 
log Q > 3, we calculate the gas opacity by extrapolating the 
tabular data. Even if the values at log Q = 3 are used for 
log Q > 3 without extrapolation, no significant difference in 
the critical core mass is found. 

The Rosseland mean opacity of dust grains is calculated 
from the monochromatic opac ities of dust grains in pro- 
toplanetary discs presented by ISemenov et al.1 (|2003h . The 
dust constituents, their evaporation temper atures, and the 
size d istribution of the dust grains used in Semenov et al.l 



Pollack et al.l 



(Il994h. 



Stognienkd 



l|2003h are the same as those in 

while their optical constants refer to lHenning 
l|l996l) . The dif ference between the grain opacities is nev- 
ertheless small jSemenov et al.ll2003l l. We introduce a fac- 
tor, /, that represents reduction or enhancement of the 
grain opacity, namely, grain opacities being / times opac- 
ities of grains in disc gas. The reasons why we intro- 
duce the factor are as follows. When icy planetesimals go 
through the upper envelope, if undifferentiated like comet 
Halley (e.g. iMumma. Weissman. fc SternI 1 19931 ) , they re- 
lease dust grains embedded in a matrix of ice upon evap- 
oration. The deposited grains raise the grain opacities in 
the upper envelope. On the other hand, opacities of dust 
grains in the envelope may be reduced due to th eir coag- 
ulation and settling as described in Introduction (IPodolal 
I2OO3.; Movs hovitz fc Pod olak 2008; Movshovi tz ct al.„201° ^ 
In this study, the total opacity is expressed by 



lay 

f 



(14) 
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Figure 2. The critical core mass (J\/crit) as a function of the mass 
fraction of heavy elements in the lower layer of the envelope, Zh- 
The solid, dashed, dotted, and dot-dashed curves correspond to 
the results for Th = 300, 500, 1000, and 2000K, respectively. 
The double dot-dashed curve is the result of wholly-polluted en- 
velopes. In all the calculations, we assume L = 1 X 10^^ erg/s and 
/ = 1. A sudden change in the slope at Zh = 0.815 for Th = 300K 
and the wholly-polluted envelopes appears because one of the ma- 
jor molecules, H2, is replaced with CO2 (see the text). 



Table 1. Free Parameters and Their Values. 



Parameter 



Value 



homopause temperature, Th 
mass fraction of heavy elements 
in the lower layer, Zh 
grain-depletion factor, / 
luminosity, L 



300, 500, 1000, 2000K 



0.015 - 0.94 
0,0.01,1,10 
1 X 10^6, 10^7, 



10^8, lO^^erg/s 



where Kgas is the Rosseland mean opacity of gas and Kg 
that of dust grains in protoplanetary discs. 



3 CRITICAL CORE MASSES 



Critic al core masses are found in the same way as iMizund 
(|l980l ). For a set of four parameters listed in Table [TJ we de- 
termine the static structure of the envelope (i.e., dS/dt = 0) 
to find the core mass for a given protoplanetary total mass, 
Afp. We increase Mp and repeat the same procedure until 
the core mass reaches a first maximum, which is the criti- 
cal core mass. In this section, we demonstrate how envelope 
pollution affects the critical core mass (Merit). Our results 
show that the heavy-element enrichment, in general, lowers 
Afcrit and this behaviour of Afcrit holds good for any choice 
of four parameters. 



3.1 Effects of Envelope Pollution by Icy 
Planetesimals 

Figure [2] shows Afcrit as a function of the mass fraction of 
heavy elements in the lower envelope, Zh, for four different 
homopause temperatures, Th. The solid, dashed, dotted, and 
dot-dashed curves correspond to Th = 300, 500, 1000, and 
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Figure 3. The critical core mass (M„it) as a function of Zj, for four different Th: (a) = 2000K, (b) Th = lOOOK, (c) = 500K, (d) 
Th = 300K. The sohd curves are the same results as Fig. 2. The dashed and dot-dashed ones exclude the rise in Kgas and the reduction 
of Vad J respectively, and the dotted one includes only the increase in 
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Figure 4. The adiabatic temperature gradient (Vad) as a function of temperature when Th = 500K and Th = 2000K. Convective and 
radiative regions are shown by solid and dashed lines, respectively. The left panel plots the results for Zh = 0.1 and the right panel plots 
those for Zh = 0.8 at M^rit. 
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2000K, respectively. We also examine the wholly pol- 
luted envelope as an extreme case, which is shown by the 
double dot-dashed curve. In all the calculations shown in 
Figure[2j we assume L = 1 x 10^^ erg/s and / = 1. 

As found in Fig. (2] except for small increases for ^ 
0.1, the heavy-element enrichment lowers Alait for any Th. 
In particular, the reduction in Merit is significant for high 
Zh. In Fig[2] one notices that Merit depends on and Th 
in complicated manners. There are three factors to change 
Merit as described in Introduction. To see how each effect 
contributes to change Merit , we have done the following sen- 
sitivity tests; the results are presented in Fig|3] The solid 
curves correspond to the results shown in Fig[2] Other three 
curves represent the results for cases where one or two out 
of the three factors are artificially excluded. 

• The dashed curves: We have used Kgas with the solar 
abundances throughout the envelope instead of including 
the rise in Kgas- The difference between the solid and the 
dashed curves represents the increment of A^ferit caused by 
the rise in fCgas. 

• The dot-dashed curves: We have used Vad with the 
solar abundances throughout the envelope to exclude the 
effect of reduction of Vad • The difference between the solid 
and the dot-dashed curves represents the decrement of 
Merit due to reduction of Vad- 

• The dotted curves: We have used both Kgas and Vad 
with the solar abundances throughout the envelope to ex- 
tract only the effect of the increase in /i. The difference be- 
tween the values of Merit for Zh ~ Za{~ 0.015) and the 
dotted curves shows the decrement of Merit due to the in- 
crease in /i. 

First, in the case of Th — 2000K, the increase in fi is the 
most important for lowering Merit, as shown in Fig|3^. The 
lower envelope (T ^ 2000K) is fully convective (see Fig|4]), so 
that the rise in Kgas has no influence on the change in Aferit . 
That is why the heavy-element enrichment in the region of 
T ^ 2000K always decreases Merit with increasing Zh- As 
for all the other cases, Merit increases at first and decreases 
beyond Zh = 0.1. When Zh < 0.1, the increase in fj, and 
reduction of Vad are inefficient in lowering Merit as shown 
in Fig[3)3-[3Ji, so that the rise in Kgas raises Merit. 

For low Th and relatively high Zh, reduction of Vad is 
the most effective in lowering Merit, as shown in Fig[3ji. To 
understand this behaviour, we present profiles of Vad in the 
envelope when Meore = Merit for Th = 500K and 2000K 
in FigH Zh = 0.1 (left panel) and Zh = 0.8 (right panel). 
The other Th cases are not presented here because the re- 
sults for Th = 300K and lOOOK are quite similar to those for 
Th = 500K and 2000K, respectively. Convective and radia- 
tive regions are indicated by the solid and the dashed lines in 
FigSl respectively. Two deep valleys of Vad are seen around 
500 - 900K and above 2000K in FigO Especially, Vad for 
500-900K is much smaller in the case of Th = 500K than in 
the case of Th = 2000K. Such small Vad is responsible for 
the small value of Merit. 

The reason for the decrease in Vad is that several chem- 
ical reactions occur in that temperature range. We present 
the number density of each chemical species, n(X), relative 
to that of molecular hydrogen, n(Il2), for two typical cases 
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Figure 6. Critical core masses (Merit) for four different / in the 
case of Th = 1000 K. Tlie solid, dashed, dot-dashed, and dotted 
curves correspond to / = 1, 10, 0.01, (no grain). In all of the 
calculations, L = 1 X 10'^^ erg/s is assume d. Only for the results 
of / = are given bv lHori fc Ikomal 1I2OIOI) . The solid curve is the 
same as the result represented by the dotted one in Fig. 2. 

in FiglSl Zh = 0.6 and Th = lOOOK (left panel) and Zh = 0.9 
and Th = 300K (right panel). Note that a sudden change in 
the slope at Zh — 0.815 of the solid line in Fig. |311 appears 
because one of major molecules, H2, is replaced with CO2 at 
Zh ^ 0.815 as shown in FigEjD. In the regime of 400 - 900K, 
there occur the abrupt decreases in H2O and CII4, which are 
stable chemical species at low temperatures, and the sharp 
increases in CO and CO2. Those chemical reactions have 
great impacts on reduction of Vad- Above 2000K, dissoci- 
ations of several kinds of molecules such as H2, H2O, and 
CO2 take place and reduce Vad- As a result, the envelope is 
rich in CO, H, and O above 2000K. 

3.2 Dependence on grain opacity 

We considered only the case of / = 1 in Section 3.1. The 
upper layer of the envelope may, however, be grain-rich or 
grain-poor, as described in Section 2. We consider both 
grain-rich (/ > 1) and grain-poor (/ < 1) envelopes in 
this subsection. The upper layer (i.e., T < Th) is convec- 
tively stable in most cases. Thus, it is to be verified that 
the spatial distribution of added dust grains in that layer is 
homogeneous; nevertheless, we adopt a constant value of / 
throughout the upper layer for simplicity. 

The behaviours of Merit for Th = lOOOK with respect 
to four different / are shown in Fig[Sl The solid curve is the 
same as that in Fig[2]D. The dashed, dot-dashed, and dotted 
curves correspond to / = 10, 0.01, and 0, respectively. Re- 
sults for / > 10 hardly differ from that for / — 10. In all the 
calculations, L = 1 x 10^^ erg/s is assumed. 

As shown in Fig. (6] Merit increases as a whole (i.e., for 
any Zh), as / increases. This is because dust grains are the 
dominant opacity source in the upper layer. The increase in 
Aferit is found to level off around / = 1. This is because 
the part of the envelope where the grain opacity dominates 
becomes convective; on the other hand, the slight increase 
for Zh < 0.1 is due to increase in gas opacity in a relatively 
high-temperature, radiative region in the layer where the gas 
opacity dominates. We have also investigated the behaviour 
of A/erit for different Th and confirmed that both cases of 



Envelope Pollution of Gas Giants by Icy Planetesimals 7 



Table 2. Critical core masses for Tjj 
four different L (erg/s) . 
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grain-rich and grain-poor envelopes show trends similar to 
that for Th = lOOOK. 
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3.3 Dependence on luminosity 

We assumed L = 1 x 10^^ erg/s in all the calculations 
shown above. The dependence of M^rit on L is qualitatively 
the same as that in the case of so lar-composition envelopes 
l|lkoma. Nakazawa. fc Emoril200oh . Here we quantify the de- 
pendence for polluted envelopes. We calculate Merit for three 
different values of L = 1 x 10^®, 1 x 10^*, and 1 x 10^^ erg/s. 
Those values are chosen based on accreti on rates of planetes- 
imals that formation theory predic ts () Pollack et al.l 1 19961 : 
iFortier. Benvenuto. fc Bruninill2009l ). The results are listed 
in Table [3 where Th = lOOOK. Table [2] shows that the de- 
creasing trend of Afcrit with is the same, irrespective of 
L. We also have found that this is true for other Th. 



3.4 Dependence on semimajor axis 

We discussed Merit at the present location of Jupiter in the 
previous subsections. In this subsection, we consider Merit 
at two different semimajor axes, lAU and lOAU. The choice 
of the two values is based on recent theories of the snow line 
in protoplanetary discs. Many relevant st udies suggest that 
the snow line moves and reaches ~ lAU ( Sasselov fc LecM 
I2OOOI : iDavisI I2OO5I : iGaraud and LinI l2007l : Win et al.ll2011^ . 
We thus adopt lAU as the innermost orbit of the protoplanet 
that can experience envelope pollution by icy planetesimals. 
The disc density (pdisc) and temperature (Tdise), namely, 
the outer boundary conditions, are taken as mentioned in 
Section 2. 

The sensitivity of Merit to the outer boundary con- 
ditions was already investigated in detail jMizuno 1980l: 
Stevensonll l"9p: IW'uchter ill 19931 : 1 Ikoma. Emori. fc Nakazawa! 



200ll : lRafikovll2006l '): Itls known that Merit for wholly radia- 



tive envelopes is almost independent of the outer bound- 
ary conditions, while Merit for fully convective envelopes 
depends significantly on them. Polluted envelopes tend to 
be convective, especially near the central star. The criti- 
cal disc density above which th e envelope is wholly convec- 
tive is given by equation (13) of llkoma. Emori. fc Nakazawa! 
(|200lf ): 



crit 
PdlBC 

Pmmsn 



VIAU/ 



3/2 



L 
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Tmmsn 
1 



5/2 



Me 



1 X 10^'' erg/s 
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V 1500K 
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Figure 7. Critical core masses (Merit) as a function of Zjr for 
Th = 300K, lOOOK, and 2000K. The same as Fig.2, but for differ- 
ent semimajor axes. The solid, dashed, and dotted curves adopt 
the outer boundary conditions of 5.2AU, lOAU, and lAU, respec- 
tively. 



where Tcvap is the evaporation temperature of grains, ft;cvap 
is the grain opacity near T = Tevap, and the subscript 
"MMSN" means that the quan tity takes the value of the 
MMSN model (jHavashil Il98l! ). Based on the fact that 

3 — 5 cm^ g~^ for other Th at Tevap 
satisfied for all the models of lAU and Th = 300K. In fact, 
the outer envelope becomes fully convective in the cases of 



i-evap 30 — 100 cm^ g ^ for Th — 300K and ftovap ^ 

lOQOK, pdise > p'St is 



(15) 



Th = 300K and ap = lAU. It is, thus, worth knowing the 
impact of outer boundary conditions on A/erit in the case of 
polluted envelopes. 

Figure [7] shows Merit as a function of at 5.2AU (solid 
line), lAU (dotted one) and lOAU (dashed one) for Th = 
300, 1000, and 2000K; the values of the other parameters 
are the same as in FiglH The trend of Merit is the same for 
all the cases. The impact of outer boundary conditions is 
small at Op > 1 AU; Merit differs at most by a factor of 2 
even for Th = 300 K and Zh ~ 1. 



4 GAS ACCRETION TIME-SCALE 

As mentioned in Introduction, gas accretion is slower for 
smaller core mass, if the other parameters being the same. 
We need to check the time-scale for gas accretion in addition 
to the critical core mass. In this section, we demonstrate that 
envelope pollution by planetesimals can also hasten the gas 
accretion on to a protoplanet. We adopt the growth time 
of the envelope, Tgas, defined in llkoma. Nakazawa. fc Emoril 
(|2000f ) as the typica l time-scale of gas accr e tion on to the 
protoplanet (see also llkoma fc Genda! (|2006h : lHori fc Ikoma! 
(2010)). We consider the protoplanet with a given core mass. 
We simulate accumulation of disc gas after the core growth 
stops. Both Zh and Th in those simulations are assumed to be 
constant in time. We discuss the validity of this assumption 
in Section 5. 

Figure [8] shows Tgas as a function of the core mass for 
Th = lOOOK, / = 1, and Zh = 0.1, 0.3, 0.6, and 0.8 (solid 
lines). For comparison, we plot the results for envelopes with 
the solar abundances: / = 0.01 (dashed line) and grain-free 
(/ = 0) envelopes (dotted one). This figure demonstrates 



8 Y. Hori and M. Ikoma 




10 



Figure 8. The growth time-scale of the envelope (rgas) as a func- 
tion of core masses. Four solid lines represent results of four dif- 
ferent Zh for Th = lOOOK and / = 1: = 0.1, Zh = 0.3, 
Zjj = 0.6, and Z^ = 0.8. For comparison, we plot the results of 
/ = 0.01 with solar abundances (dashed line) and grain-free en- 
velopes with the solar abundan ces, / = (dotted one), which are 
found bv lHori fc Ikomal | |2010|V 
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Figure 9. Core masses with Tgas = IMyr for three : Tji = 
lOOOK (sohd), 500K (dashed), 300K (dotted). In all the calcula- 
tions, we assume / = 1. 



Table 3. Core masses with Tg; 



IMyr. 



that envelope pollution by icy planetesimals accelerates gas 
accretion on to the protoplanet significantly. The accelera- 
tion is more effective compared to that by reduction of grain 
opacities in the envelope. We do not suppose reduction of 
grain opacities in the upper envelope in those simulations. In 
addition to envelope pollution, if we take into account reduc- 
tion of grain opacities, it gives positive feedback to shorten 

gas 

(see Table El). 

Gas giant formation must be finished before disc gas 
disappears. The disc lifetime is believed to be 1-10 Myrs 
based on the disc frequency estimated from near-IR ex- 
cess (JHKL-hand) observations of stars in young clus - 
ters (jHaisch. Lada. fc Ladal I2OOII : iHernandez et al.l I2OO7I ). 
We now consider the core masses with rgas ^ IMyr as allow- 
able for gas giant formation. Minimum core masses (M™!.^) 
that satisfy the criterion (rgas ~ IMyr) are listed in Table 
[3] for f — 1 and / = 0.01. We cannot simulate quasi-static 
evolution of the protoplanet 's envelope with smaller cores 
than O.IM0 because of numerical difficulty. In such a case, 
we show M™j° < O.IM® in Table [3 In FigH we also plot 
M™ as a function of for three Th, lOOOK (solid), 500K 
(dashed), and 300K (dotted). The behaviours of ACrS with 
respect of are quite similar to those of Merit- We find that 
envelope pollution has contributions to hasten Tgas in the 
similar manner to the lowering of Ma-it- Our results suggest 
that envelope pollution by icy planetesimals has the poten- 
tial to make gas giant formation with small cores possible; 
for example, gas giants with cores smaller than IMg can 
capture disc gas by IMyr when Zh ^0.7 for / = 1 or when 
^ 0.5 for / = 0.01. 



5 DISCUSSIONS 

5.1 Dilution during runaway gas accretion 

We assumed that both Zh and Tk are constant with time. 
This may be oversimplification and questionable, especially 
in the phase of runaway gas accretion. When the critical core 
mass is attained, the accretion rate of disc gas is much higher 
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than that of planetesimals. In an extreme case where the un- 
polluted outer envelope never exchanges material with the 
polluted lower envelope, the gas accretion results in increas- 
ing the mass only of the upper envelope. The lower envelope 
behaves like a part of the "core" . In this case, the envelope 
pollution does not resolve the problem of the slow formation 
of gas giants with small cores. 

In the other extreme case where the inner and outer 
envelopes exchange material instantaneously between each 
other via eddy diffusion, accreting fresh disc gas dilutes the 
polluted lower envelope, which results in decelerating the 
disc-gas accretion. On the other hand, the dilution is in- 
evitably accompanied by mass growth of the envelope, which 
accelerates the disc-gas accretion. 

Unfortunately, it is uncertain whether or not mixing 
occurs effectively, as follows. A characteristic time-scale of 
eddy diffusion, Teddy, is given by 

T-cddy ~ Jl^, (16) 
22 

where Hp is the pressure scale-height and 7^22 is the coeffi- 
cient of eddy diffusion. Estimated values of Teddy are listed 
in Table 4, where we have used our numerical values of Hp 
at the tropopause when Mcore = Merit. We have adopted 
K,, = 10^ - 10** cm s which c orrespond to values for 
the p resent Jupiter and Saturn (e.g. lChamberlain fc HuntenI 
[l983), although K^^ for protoplanetary envelopes is uncer- 
tain. The lower and upper values of Teddy at each Z-t and K^z 
in Table 4 are Teddy of Th = 2000K and that of Th = 300K, 
respectively. As seen in this table, Teddy ranges widely from 
10^ to 10* yr. Since time-scales of gas accretion of interest 
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Table 4. r. 
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at the critical core masses. 
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are 10^-10^ years (which correspond to the observed disc 
lifetime), we can derive no definite conclusion as to whether 
dilution occurs or not, because of the uncertainty of Toddy 
Mixing in accreting protoplanetary envelopes is found to be 
an important issue with gas giant formation. 



5.2 Non-ideality of envelope gas 

In this study, we considered the envelope gas as ideal. In the 
case of the envelope with the solar abundances, we find that 
non-ideal effects cause less than 10% change in the critical 
core mass. This paper aims to point out the importance 
to incorporate the effects of heavy-element enrichments in 
simulations of the structure and evolution of protoplanetary 
envelopes and to sort out what has a great impact on the 
critical core mass. Including the effects of non-ideality of gas 
is a future study. We found that reactions between molecular 
compounds play a important role in determining the critical 
core mass. Thus, we need thermodynamic quantities of non- 
ideal mixture of at least H, He, C, and O. 



5.3 Dynamical stability 

We should al so pay attention to the dyn amical stability of 
the envelope. iTaiima &c Nakagawal (| 19971 ) found that the en- 
velope with solar abundances is stable dynamically through- 
out quasi-static evolution up to Jupiter's mass. However, it 
is not obvious whether the envelope polluted by heavy ele- 
ments is stable or not. In this study, reduction of Vad due to 
molecular dissociations is one of important effects to lower 
Merit and shorten Tgas . Vibrational instability due to molecu- 
lar dissociation may excite nonlinear hydrodynamical waves 
by K- mechanism bec ause of reduction of Vad and the rise in 
Kgas (|Wuchter]||l99ll ) . Nonlinear hydrodynamical waves may 
cause dynamical expansion of the envelope and the ejection 
of a large part of the envelope. Thus, the dynamical stability 
of the polluted envelope should be also investigated in the 
future. 



5.4 Implication for Jupiter and Saturn 

It is worth while to consider the origin of Jupiter and 
Saturn from the viewpoint of core masses. Interior mod- 
ellings of the two planets suggest that Jupiter's core is 
smaller than lOMffi, while Sat u rn's core is larger tha n 
lOMffi dg aumon fc Guiiio^ |2004| : iNettelmann et al.l l2008t ). 
iMilitzer et al.. (.20081 ') reported that Jupiter may possess 
a massive core of > lOM®, but their interior modellings 
with rigid rotation are inconsistent with the observed value 
of the 4th-order gravi t ationa l moment, J4. In addition, 
iFortnev fc NettelmannI (|2010l ') pointed out that the mass 



fraction of helium used by IMilitzer et al] (|2008l l is responsi- 
ble for a massive core of Jupiter. Although the exact mass 
of Jupiter's core is still an open question, the core mass 
of Jupiter suggested by interior modellings is, on average, 
smaller than that of Saturn. 

In this study, we demonstrated that envelope pollution 
by icy planetesimals enables formation of gas giants with 
small cores, which seems compatible with recent modellings 
of Jupiter. If so, how about the origin of Saturn that is in- 
ferred to have a large core? The large core inside Saturn may 
reflect the history of high accretion rates of planetesimals. 
A variety of ideas for high accretion rates have been inves- 
tigated so far — e.g., the in fluence of growing proto- Jupiter 
and t he other protoplanets (|Guilera. Brunini. fc Benvenutd 
201(]|) . Type I migration of protoplanets ( Alibert et al.l 
200a), hi gh stellar m etallicity and a large initial disc mass 
(Tliommos fc Murray 2006), high initial surface density 
of so lid materials (|Dodson-Robinson et al.l l2008l : iRaflkovl 
l201lf ). However, high accretion rates of planetesimals could 
have a risk of forming Saturn with a small core. This is 
because a large amount of small-sized fragments are pro- 
duced in the Saturn's region due to frequent high-velocity 
collisions of planetesimals, so that the envelope should be 
polluted heavily by icy planetesimals. Thus, we require to 
revisit accretion rates of planetesimals on to proto-Saturn, 
for example, in the case where Jupiter flrst emerges and 
other planetary cores still continue to grow. This is because 
planetary accretion proceeds more rapidly in the inner re- 
gions, so that proto- Jupiter may have formed a small core 
flrst and triggered runaway gas accretion due to envelope 
pollution by icy planetesimals. 

As one alternative idea for a large core of Saturn, 
iLi. Agnor. fc LinI (|2010l ) suggested that Saturn might have 
experienced an impact of a (proto-) gas giant and then 
merged. It is still an open question whether such an impact 
on Saturn happened in reality in a multiple-protoplanet sys- 
tem. Direct N-body simulations of competitive core growth 
beyond the snow line will give an answer to this problem. 

We should also comment on another possibility that 
the disc instability scenario may be able to explain gas gi- 
ants with small cores. The disc instability scenario is usu- 
ally a single-stage model in which a massive protoplane- 
tary disc, typically O.IMq, becomes gravitationally unstable 
and rapidly co l lapses to form a proto-gas giant planet (e.g. 
Cameronlll97i; IBoss|[200^: iMaver et al.ll2007l: iDurisen et al] 



20071 : lBolevll2009l : llnutsuka. Machida. fc Matsumotoll2010l ) 



Recent works on the disc instability scenario showed that 
a gaseous clump is able to capture planetesimals in the 
course of its rapid contraction, assuming concurrent forma- 
tion of large-sized planetesimals (iHelled. Podolak. fc Kovetj 
l2008l : iHelled fc Schubert! I2OO8I . |2009| ) or a gaseous clump 
is formed in a dense region where solid materials (1/^m- 
10cm) are collected aerody namically by spiral arms 
iBolev. Helled, fc Pavn3 I2OI1I ). In those scenarios, grains 
that grow to typically ~ 10 cm can settle down toward the 
centre against the stirring of convective motion and form 
a relatively small core. Also the disc instability scenario 
is faced with the problem of Saturn formation, although 
Jupiter with a small core could be formed by the disc in- 
stability. One possibility to form a large core of Saturn is 
that the existence of Jupiter influences the distribution of 
planetesimals near Jupiter's region and enhances the core 
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growth of Saturn. As mentioned in the above paragraph, 
this situation is the almost same as the envisioned picture 
of the Saturn formation in the core-accretion m odel taking 
account into envelope pollution. More recently, iNavakshinI 
has proposed another idea of gas giant formation 
with large cores by the disc instability: Grain sedimenta- 
tions inside a large gaseous clump yield a large core and 
trigger rapid gas accretion on to the core. When the clump 
migrates into a few AU from a parent star, it experiences 
tidal disruption and evaporation due to stellar irradiation. It 
is peeled off the outer metal-poor envelope of the clump and 
then results in a Saturn-like planet. In the future, further 
detailed studies will verify his idea. 

5.5 The amount of heavy elements in envelopes of 
extrasolar gas giants 

Recent measurements of the mass and radius of transiting 
extrasolar gas gia nts (EGG) allo w us to constrain planetary 
composition (e.g. iGuillotI 120081 ). Those internal structures 
showed a variety of the total amount of heavy elements con- 
tained in transiting EGGs, including cores. Although it is 
not clear whether such a diversity of the amount of heavy el- 
ements in EGGs reflects a diversity of their core masses, core 
masses of EGGs are a key test to verify how efflciently our 
idea of envelope pollution works in reality. Further transit 
observations of EGGs and our comprehensive understanding 
of interior structures of EGGs will help us to reveal forma- 
tion of giant planets. 

5.6 Implication for Uranus and Neptune 

Finally, we consider an implication of interior structures of 
Uranus and Neptune. Envelope pollution by icy planetesi- 
mals (and erosion of their icy shells) may be in good agree- 
ment with interior modellings of the present Uranus and 
Neptune. The H/He envelopes enriched by volatile molecules 
such as H2O and CH4 have been proposed in order to re- 
produce their observed gravi tational moments (J2 and J4), 
which are poorly constrained (iPodolak^ Weizman. fc MarlevI 



results in the polluted H/He envel ope because the current 
Uranus has the extreme axial til t (|Korvcanskv et al]|l990l : 
ISlatterv. Benz. fc Cameronlll993 ). 



ll995l : lMarlev. Gomez, fc Podolaklll995^ . Manv recent works 
on interior modellings of Uranus and Neptune also suggested 
that they may have the H/He envelopes into which "ices" 
are mixed instead of traditio nal three-layer models, that 
is density-stratified inte riors (|Podolak. Podolak. fc MarlevI 
I2OOOI : iHelled et al.|[2oTll ). On the other hand, the fact that 
the present Uranus and Neptune have non-axisymmetric, 
non-dipolar magnetic fields may favour the existence 
of stably-stratified fluid layers bene ath the outer thin 
shell that drives dynamo action (|Stanlev fc BloxhamI 
I2OO4I . I2OO6I ). In addition, the compositional stratiflcation 
in the interior m ay be re sponsib le for extremely low- 
luminous Uranus jPodolak. Hubbard, fc Stevenson Il99ll : 
iHubbard. Podolak. fc StevensonI Il995l ) It is still unclear 
whether the H/He envelopes are gradually- mixed with ices 
toward the centres or not. In any case, both Uranus and 
Neptune are likely to possess the outer polluted envelopes 
from the viewpoint of the gravitational harmonics and ther- 
mal evolution. Their polluted envelopes may be consistent 
with the picture of envelope pollution by incoming icy plan- 
etesimals, although we cannot rule out the another possibil- 
ity for Uranus that a giant impact on a primitive Uranus 



6 SUMMARY AND CONCLUSIONS 

We have investigated how envelope pollution by incoming 
icy planetesimals affects critical core masses (Merit) and 
growth time-scales of the envelope (rgas). We considered that 
the envelope has two-layer structure: the upper layer with 
the solar abundances (non-polluted layer) and the lower one 
with mixture of solar and comet-Halley-like compositions 
(polluted layer). We introduced two key parameters, the 
mass fraction of the heavy elements in the lower envelope 
(Zh) and the homopause temperature (Th) which defines 
the boundary between the two layers. The main results that 
we obtained in this study are summarized as follows. 



(1) Critical core masses: 

Envelope pollution by icy planetesimals lowers Merit 
for any choice of values of Th and most of the range of 
Zh{> 0.1) (Figig. Widely-polluted (low Th) and highly- 
polluted (high Zh) envelopes lessen Merit remarkably. The 
increases in the molecular weight (/x) and reduction of 
adiabatic temperature gradient (Vad) are responsible for 
the lowering of Merit (Fig. In particular, remarkable 
reduction of Vad due to chemical reactions for 400-900K 
and molecule dissociations for > 2000K (Figs. |4] and [5]). 
Such behaviour of Merit is the same for any choices of other 
parameters, the grain opacities (Fig. |6}, the luminosity 
(Table [2]), and the semimajor axes (Fig. [Tjl. 

(2) Growth time-scales of the envelope: 

We have simulated the growth of the protoplanetary 
envelope after the core growth stops. Our simulations 
demonstrate that the envelope pollution can also hasten 
gas accretion on to the protoplanet (Fig[8]), which implies 
that envelope pollution may enable gas giant formation 
with small cores (Fig. |9] and Table |3} . 



To confirm that our idea really works, further studies 
will be required: (i) chemical evolution in the polluted enve- 
lope coupled with eddy diffusion and infiow of disc gas, (ii) 
dynamical stability of the polluted envelope, (iii) competi- 
tive core growth (planetary accretion) beyond the snow line 
in a multiple-protoplanet system, (iv) the efficiency in mass 
deposits of icy planetesimals based on their trajectories and 
sizes in the envelope, (v) envelope pollution due to erosion of 
a core itself, and (vi) chemical equilibria of non-ideal gases 
in the polluted envelope. However, it is quite certain that 
envelope pollution by planetesimals occurs and affects the 
formation process of gas giants. We claim that it is necessary 
to take into account envelope pollution when we discuss gas 
giant formation based on the-core accretion model. 
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